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Abstract 

We have studied the electron-phonon and superconducting properties of the Mgi_2;Alx-B2 and 
^^2{i-y)^2y alloys within the framework of density functional theory using the self-consistent 
virtual-crystal approximation. For both alloys, the Eliashberg spectral functions and the electron- 
phonon coupling constants have been calculated in the two-band model for several concentrations 
up to x{A\)= 0.55 and y{C)= 0.175. We solved numerically the two-band Eliashberg gap equations 
without considering interband scattering. Using a single parameter for the Coulomb pseudopoten- 
tial, which was determined for the undoped compound, we were able to reproduce the experimental 
doping dependence of A^, A^r, and Tc for both alloys on a quantitative level. In particular, the 
observed differences in the doping range of superconductivity between Al and C doping indicate a 
pronounced influence of the doping site, which can be explained naturally in the present approach 
without the need to invoke interband scattering, suggesting that this factor plays only a minor 
role. 

PACS numbers: 74.70.Ad, 74.62.Dh, 63.20.kd, 71.15.Mb 
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I. INTRODUCTION 



The discovery of superconductivity in 2001 in the intermetalhc compound MgB2, with 
a critical temperature Tc ~ 39 has motivated a lot of theoretical and experimental 
studies in order to understand the origin and characteristics of the relatively high Tc in this 
material. It is now generally accepted that MgB2 is a phonon-mediated superconductor, 
and that the high transition temperature arises due to a combination of several peculiar 
features in its electronic structure and electron-phonon coupling, which conspire to produce 
a superconducting state with multiple gaps^i"— . Its electronic band structure in the vicinity 
of the Fermi energy consists of two bonding cr bands corresponding to in-plane s — Px — 
Py (sp^) hybridization in the boron layer and two vr bands (bonding and anti-bonding) 
formed by hybridized boron pz orbitals. A substantial part of the electron-phonon coupling 
has its origin in the interaction of states at the cr-band Fermi surfaces with one specific 
phonon mode, the B-B bond stretching mode with E2g symmetry at the F point^"^*^"— . In 
addition, MgB2 possesses two distinct superconducting gaps associated with the a and vr 
Fermi surfaces. This superconducting state can be described within a multiband version 
of the Eliashberg theory where the pairing interaction is split into intra- and interband 
contributions^iii^"— . 

In the search for related compounds with similar outstanding superconducting properties, 
only a few variants of MgB2 have been found. Among them, two alloy systems could be 
successfully synthesized based on the partial substitution of Mg by AU^"— and B by (J^-^^ 
respectively. Both substitutions provide electron doping to the alloy and lead to a reduction 
of Tc. For Mgi_2;Alx-B2 and MgB2(i_j,)C2y, loss of superconductivity is found for x > 0.5^"— 
and y > 0.15^"—, respectively. A correlation between the Tc reduction and the filling of the 
hole-type cx-bands as a function of Al-doped content was found by first principles calculations 
within the virtual- crystal approximation2Ii2^. In parallel with the reduction of Tc also a 
decrease of the superconducting gaps Ao- and A^r has been observed. For single crystals as 
well as polycrystals Al doping was found to decrease both a and vr gaps monotonically, which, 
however, stay distinguishable even for Tc as low as 10 K (x ~ 0.32)^"—. These observations 
indicate that the interband scattering F^^ remains small even at high doping levels, and 
is insufficient to produce a merging of the a and vr gaps. On contrast, for the C-doped 
system, contradictory experimental results have been reported with respect to the question 
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if the superconducting gaps merge as a function of doping^i^i^"— . In all experiments a 
decrease of both gaps with C doping was observed. There are point-contact tunneling^^i^, 
point-contact spectroscopy^^ and photoemission spectroscopj*^^ measurements that show a 
clear difference between the a and vr gaps at all doping levels. However, there exists also 
point-contact spectroscopy measurement s^i that suggest a merging of the gaps at Tc ~ 17 
K 0.13). This was then interpreted as a doping-induced increase of the interband 

scattering Fo-tt, which tends to reduce gap anisotropics. 

From the theoretical point of view, many investigations have been performed to study 
the doping dependence of the structural2Ii2£^ electronio^^'^"— , vibrational^"— and supercon- 
ducting properties^'^'^'^"— of the Mg1_a.Ala.B2 and MgB2(i_j;)C2y systems using different 
approximations for the simulation of the alloys, like the supercell approacb^^'^'^'^'^, the 
rigid band approximation (RBA)^, the virtual- crystal approximation (VCA)^i^i^i^, the 
coherent potential approximation (CPA)^i^i^, and the Korringa-Kohn-Rostoker coherent 
potential approximation (KKR-CPA)^. However, in particular for the supercell and CPA 
approaches, the studies have been limited to a few Al or C concentrations only, because 
these calculations are computationally very demanding, especially if one is interested in 
very low (close to Mg or B) or high (close to Al) concentrations. Additionally, in these ap- 
proaches the symmetry of the original system is lost, which complicates the interpretation 
and understanding of experimental results as a function of doping. 

In this paper, we present a study of the electron-phonon and superconducting properties 
of Mgi_a:Ala.B2 and MgB2(i_2/)C2y within the framework of density functional theory^ using 
the self-consistent virtual-crystal approximation (VCA)^i^"— . We calculate the electron- 
phonon (e-ph) properties such as the Eliashberg function, afjF{uj)^ and the e-ph coupling 
constant, Ajj, within the two-band model as a function of doping. By solving the two-band 
Eliashberg gap equations on the imaginary axis we obtain the superconducting gaps, A^. and 
Att, and the value of as a function oi x 01 y for the Mg1_a.Ala.B2 and MgB2(i_j,)C2y alloys, 
respectively. The evolution of these quantities is analyzed and discussed in connection with 
changes in the electronic and vibrational properties. 
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II. COMPUTATIONAL DETAILS 



The calculations were performed with the mixed-basis pseudopotential method (MBPP)^'^. 
For Mg/Al and B/C norm-conserving pseudopotentials were constructed according to the 
Vanderbilt description^. Details of pseudopotentials, basis functions and calculational as- 
pects for ground-state and phonon properties can be found in a previous publication^^. 
The Mgi_a;Ala;B2 and MgB2{i_j;)C2y alloys were modeled in the self-consistent virtual- 
crystal approximation (VCA)^'2^'^"— i^"— . The VGA is implemented within the MBPP 
method^i^ by generating new pseudopotentials with a fractional nuclear charge at the 
Mg or B site for each x and y, respectively (Al: Z=12+x and C: Z=5+y), and by ad- 
justing the valence charge accordingly^^. From our previous results for the electronic and 
vibrational properties^ the screened electron-phonon matrix elements were calculated via 
density functional perturbation theory^i^"— , which are the key elements of the Eliashberg 
theory^"—. The calculations employ the PBE version of the generalized gradient approx- 
imation (GGA)^"— , and are performed at the optimal lattice parameters for each doping 
level^. Eliashberg functions for all band combinations where obtained by standard Fourier 
interpolation of quantities calculated with a dense 36 x 36 x 36 /c-point mesh and a 6 x 6 x 6 
g-point mesh. The original four-band Eliashberg functions are projected onto an effective 
two-band model by averaging over the two a and the two vr bands, respectively. The partial 
and total Eliashberg functions are given by the following expressions, 

* qi/ k,k„ 

where i and j are the band indices a or tt, Ni (Ntot) is the partial (total) electronic density 
of states at the Fermi level (per atom and spin), and 5'k'iknj e-ph matrix element 

for scattering of an electron from a Bloch state with momentum k to another Bloch state 
k„ = k + q by a phonon qy {y indicates the branch index and Wqj, is the phonon frequency). 

In a similar way, the partial and total e-ph coupling parameters (A) are expressed as 
follows, 
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A,, =2 j —alF{u), 



(3) 



Hot ^ 



Within the two-band model, there are three independent contributions to two 



intraband (tttt and aa) and one interband al^F{uj) = jf-al^F{uj). With the knowledge of 
Q!? F(u;), the two-band Eliashberg gap equations^"— >^ on the imaginary axis were numer- 
ically solved in order to obtain the gap values and for each given Al or C concentra- 
tion, respectively. This procedure has been previously used in similar studies of undoped 
MgB^^i^. The solution involves four non-linear coupled equations for the Matsubara gaps 
Ai{iUn) and the renormalization functions Zi{iUn), 

Ai{iUn)Zi{iuJri) = TtT E [Aij{iuj,rn - i^n) " IJ^*j{uJc)0{uJc " 1^^!)] N]^i{iUJm) , (5) 

Zi{iun) = IH Aij{ium - iUn)Nj^Q{iu^), (6) 

— 

where 6 is the Heaviside function, and fi*j is the Coulomb pseudopotential, Uc a cutoff 
frequency (chosen as Uc ~ lOw^"^) and Un = vrT(2n — 1), with n = 0, ±1, ±2, . . ., is the 
discrete set of Matsubara frequencies. The pairing interaction is contained in the kernel 

Aij{lU]m -tUJn) = 2 / ; r^, (7) 



and we defined the following quantities: 



m 



(8) 



Nioi^ur^) = I (9) 

In order to keep the number of adjustable parameter to a minimum, we approximated 
the Coulomb pseudopotential matrix which is a two-by-two matrix in the case of the 
two-band model, by a diagonal form proposed earlier— as yU*^ = fi^Sij.— The gap values 
were identified with Ai{iui), which corresponds to the point on the imaginary axis which is 



closest to the real axis. Test calculations solving the Eliashberg equations on the real axis 
indicated that this approximation is accurate on the level of 1% or better. 



III. RESULTS AND DISCUSSION 



Based on our previous results for the electronic and vibrational properties=^^ we have 



calculated the electron-phonon coupling quantities of the two-band model {afjF{u) and A 



with i,j = cr, vr) for the ranges x ^ 0.55 and y ^ 0.175 in the Al- and C-doped systems, 
respectively. 
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FIG. 1. (Color online) Evolution of the (a) total Eliashberg function and (b) components a'fjF{uj) 
= (T,7r) for the Mgi_3;Al3;B2 alloy. 

In Fig. 1 we show the evolution of the Eliashberg functions, the total spectra and the four 
components afjF{u) = cr, vr) of Mgi_j;Alj;B2 for six Al concentrations in the supercon- 
ducting regime (x =0.0, 0.1, 0.2, 0.3, 0.4, and 0.55). We observe that for MgB2 (x = 0) the 
largest contribution to the total spectral function comes from the aa component, where the 
main peak centered at approximately 70 meV corresponds to frequency of the £^2g-phonon 
mode. The vrvr spectrum has its main contribution from the high-frequency phonon region 
related to the i?ic,-phonon mode, while the interband contribution, an{7ia), is concentrated 
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in the region between 50 and 70 meV. We note that, although the aa part represents the 
main contribution to the total spectra, the other components can not be neglected in a 
proper quantitative description of the e-ph coupling and of the superconducting properties. 
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FIG. 2. (Color online) Evolution of the (a) total Eliashberg function and (b) components afjF{u}) 
= o-,7r) for the MgB2(i_y)C2.y alloy. 

From the evolution of spectral functions for Mgi_^Ala;B2 (Fig. 1) we observe that almost 
all components are reduced by Al-doping, but the largest changes are exhibited by al^F{u). 
Its main peak shifts to higher frequencies and its area is reduced at the same time until it 
practically vanishes for x = 0.55. This doping level is close to the region where the loss 
of superconductivity has been observed experimentally {x > 0.5)i^"— . The reduction of 
a'^^F{uj) indicates the loss of intraband e-ph coupling between the a states and the bond- 
stretching phonon modes and has its origin in the continuous filling of the a bands, which is 
completed at Xc = 0.57—. The shift to higher frequencies is due to the hardening of the i?2g- 
phonon mode as x increases, a phenomenon discussed previously^. Similar to aa, the na 
interband contribution is also reduced as a function of x and almost disappears at x = 0.55. 
On contrast, the an contribution shows an slight increase around 50 and 70 meV, and the 
nn contribution at higher frequencies strengthens slightly with doping, but the position of 
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its main peak is almost unaffected. In recent electron tunneling spectroscopy measurements 
on Al-doped thin films^^, this general behavior of the Eliashberg function indeed has been 
observed, supporting our results. 
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FIG. 3. (Color online) Evolution of \ij and Ni{Ep) as a function of x and y for Mgi_a;Al2.B2 and 
MgB2(i„y)C2y, respectively. 



In Fig. 2 we present the results for the Eliashberg functions of MgB2(i-j,)C2j/ for six C 
concentrations in the superconducting region of the alloy {y =0.0, 0.05, 0.1, 0.125, 0.15, 
and 0.175). The different components of a^F{u!) exhibit the same trends with increasing 
C concentration as those found for the Al-doping. However the shape of the spectras are 
different and the changes take place at lower concentrations. When comparing the two 
alloys, one should take into account that the number of doping-induced electrons per unit 
cell is given by x and 2y, respectively. Even with this factor of two, the vanishing of the 
(Tcr and an contributions at 2y ^ 0.35 occurs at a much smaller doping level than for Al- 
doping (x = 0.55). The dramatic reduction of a'^^F{u!) at ?/ ~ 0.175 correlates also with 
the complete filling of the a-bands on MgB2(i_y)C2y— . 

In Fig. 3 calculated total and partial contributions for A as well as for N{Ep) are shown. 
For MgB2 the calculated values are Ag-o- = 0.850, Aq-tt = 0.196, Xt.^ = 0.145, Aj^.^- = 0.250, and 
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Atot = 0.672. It is worth mentioning that this Xtot value is very close to the experimental 
one by Geerk et ai— [Xeff = 0.650). N{Ef) partial contributions for MgB2 are N^r = 0.148 
states eV~ Vspiii and A''^ = 0.200 states eV~ Vspin, which are very similar to those calculated 
earlier by Liu et ai— and Golubov et al.—. As seen from Fig. 3, the main contribution to 
the e-ph coupling (A) in undoped MgB2 comes from the aa component. Among the different 
contributions of the e-ph coupling, Ao-o- shows the largest changes on doping with a reduction 
of ^ 75% (comparing the boundary concentrations). The other components also decrease 
with doping, albeit at different scales, ranging from the nearly constant behavior of Ao-tt and 
Atttt to an almost complete vanishing of A^^-- As a consequence, Xtot monotonically decreases 
with doping. 
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FIG. 4. (Color online) Evolution of Xij/Nj and Xtot/^tot as a function of x and y for Mgi_a;Ala;B2 
and MgB2(i_y)C2j/, respectively. Note that XT^a/^a = Xan/^^-K- 



Doping-induced changes in the coupling constants can arise from changes in the partial 
density of states or from changes in the e-ph matrix elements. In order to distinguish 
between these two possibilities, we plotted in Fig. 4 the ratios Xij/Nj and Xtot/Ntot- Indeed, 
a relationship afjF{uj) ~ Nj and Xij ~ A'^- can easily be derived from Eq. 1 under the 
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assumption of momentum- independent e-ph matrix elements. The ratios \ij/Nj for the 
interband (Ao-tt, A^o-) as well as for the intraband A,r7r couplings remain practically constant as 
a function of doping for both alloys, indicating that the corresponding e-ph matrix elements 
are approximately independent of doping. However, for the (acr) intraband coupling, the 
ratio exhibits a stronger variation with doping, in particular for Al doping, which signals a 
clear doping dependence of the e-ph matrix elements. In this case, a simple scaling with the 
partial density of states would be inappropriate to describe the doping dependence of Xaa- 




T(K) 

FIG. 5. (Color online) Temperature dependence of the superconducting gaps Aq- and (solid 
lines) for undoped MgB2 as obtained from the two-band Eliashberg gap equations. Symbols rep- 
resent experimental data (a)^, (o)^, (EB)^, and (♦)^'2i. The calculated gap values at T ^ K are 
A,^ = 7.04 meV and A^ = 2.71 meV. 

To solve the Eliashberg gap equations, we determined the single remaining parameter, the 
Coulomb pseudopotential /Xq, by the requirement that for undoped MgB2 the experimental 
transition temperature of Tc = 38.82K— is reproduced. For a cutoff frequency Uc = lOa;^"^ 
we found /Xq = 0.107. The resulting temperature dependence of the superconducting gaps 
is shown in Fig. 5 and compared with available experimental data^"-"^. The gap values for 
T — K are A^- = 7.04 meV and A^^ = 2.71 meV, respectively, in good agreement with 
experimental results^i^i^. 
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FIG. 6. (Color online) Superconducting gaps A^- and A^^ at T — t- OK and critical temperature 
Tc for (a) Mgi_a;Al2:B2 and (b) MgB2(i_y)C2j; as a function of x or y, respectively. The lines 
are the present calculations and the symbols represent various experimental measurements; (a): 
(M)^, (o)20, (*)30, (v)^, (A)^, (ffl,H)^^, (*}^; (b): (x)^, (*)^^, (.)^, (f)^^, (v)M, 



The same two-band Eliashberg approach was adopted for the alloys Mgi_2.A^B2 and 
MgB2(i-y)C2j/ keeping /Iq at the obtained value for undoped MgB2. The doping dependence 
of Act, Att, and Tc for both alloys is presented in Fig. 6 and compared with experimental 
dateiF^^^—^'—^^. The calculations reproduce the experimental trends that both gaps 
and Tc decrease with increasing Al or C doping. Beside this common feature the two alloys 
exhibit also striking differences. The first concerns the doping range where superconductivity 
exists. Tc goes to zero close to the critical concentration for which the cr-band is completely 
filled (xc(Al)= 0.57, yc{C)= 0.177)^. Thus, superconductivity vanishes significantly faster 
on C doping than on Al doping, even when taking into account that one should compare 
doping levels x = 2y as discussed above. The second difference relates to the shape of the 
Tc versus doping curves. For Al doping, Tc initially drops fast and develops a longer tail, 
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whereas for C doping Tc is only slowly reduced initially, while it exhibits a steeper drop 
towards the critical concentration where Tc vanishes. A similar difference in shape is also 
observed for the larger gap. As both alloys are electron-doped systems, these differences 
indicate the importance of the doping site for the superconducting properties. As explained 
the origin of this difference can be traced back to the distribution of 



in Refs. 



27 and 



28 



the extra charge introduced by doping. In the Al-doped system an important portion of the 
extra electrons is located in the interplanar region, and only a small fraction in the boron 
planes. In contrast, for C doping the extra charge mainly remains in the area between the 
B atoms within the boron plane, exactly in the region of the cr-bonds. Therefore, the extra 
charge introduced by C doping is more effective in reducing the number of holes in the a- 
band and has a stronger influence on the phonon frequencies, in particular, on the hardening 
of the mode. Consequently, C doping leads to a faster decrease of the e-ph coupling and 
of the superconducting properties as compared to Al doping. 

The various experimental data sets for and the gaps plotted in Fig. 6 exhibit a 
clear spread indicating a large dependence of the superconducting properties on the sample 
preparation methods and on the physical conditions of the measurement procedure itself. 
In addition, an accurate determination of the actual doping concentration in these alloys 
is complicated and far from trivial. Furthermore, there is so far no consensus about the 
behavior of the gaps for larger C doping. While some experiments suggest a merging of the 
Act and A^r gap at ?/ ~ CIS^"^, others find two distinct gaps even for the highest doping 
levels^*^"^. Within these experimental uncertainties, our calculations agree quantitatively 
with the data for both alloys. In particular, the different doping regimes are obtained in a 
natural way. We recall that our study involved only a single free parameter, /Iq, which was 
fixed for the undoped system and which does not directly affect the doping dependence or 
the gap anisotropy. 

In agreement with experimental data, the present calculation predicts for both alloys a 
stronger influence of doping on the a gap, which follows approximately the doping depen- 
dence of Tc. On contrast, the vr gap remains rather stable and only slowly decreases on 
doping. This is at variance with a previous ab initio study based on the fully anisotropic 
gap equations by Choi et al.—, where doping was modeled by simply introducing excess 
electrons. For a moderate doping level of x = 0.2 {y = 0.1), they found a severe degradation 
of the 71 gap while the a gap was more robust. This failure of a rigid-band-like type approach 
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indicates that a more self-consistent site-dependent treatment of the doping is required for 
a proper description of the superconducting properties in doped MgB2. 

Two previous computational studies^*^ of the superconducting properties of MgB2 alloys 
adopted a scaling scheme to describe the doping dependence. The Eliashberg functions for 
the undoped compound were scaled taking into account the doping dependence of N{Ep) 
and of the phonon frequency. As such an approach does not discriminate between the 
doping sites, Kortus et alr^ argued that the differences observed for Al and C doping are 
due to a larger interband scattering for C than for Al doping. The present study, however, 
demonstrates that the difference between Al and C doping appears naturally within the 
VGA approach, without the need to introduce another free parameter such as the interband 
scattering, as long as the influence of doping on the structure and on the lattice dynamics 
is properly taken into account. 

IV. CONCLUSIONS 

We have performed a first-principles study of the electron-phonon coupling and super- 
conducting properties for the Mg1_a.Al2.B2 and MgB2(i-y)C2j; alloys as a function of x and 

respectively, by combining the self-consistent virtual-crystal approximation and the two- 
band Eliashberg model. For undoped MgB2, the Eliashberg function possess a main peak 
at around 70 meV related to the -E2g-phonon mode coming from the aa contribution, and 
a sharper peak at 90 meV which originates largely from the tttt contribution, and is related 
to the i?ig-phonon mode. The total coupling constant Xtot = 0.67 agrees with the experi- 
mental value of 0.65 as deduced from tunneling measurements. For the alloys we found that 
a'^F{uj) depends very sensitively on doping. It exhibits pronounced changes both in shape 
and in position of its main peaks, which renders any attempt to derive it from the spec- 
trum of undoped MgB2 via scaling procedures very unreliable. The calculated evolution of 
the Eliashberg functions compare well with recent electron tunneling spectroscopy measure- 
ments on Al-doped thin films. The e-ph coupling parameter and its different contributions 
decrease as a function of doping for both alloys. Although both Al and C dopants donate 
electrons, the e-ph coupling exhibits a clear dependence on the doping site, which is also 
reflected in a^F(cj). With the Coulomb pseudopotential fixed for the undoped compound, 
we could reproduce the experimental doping dependence of A^., A^, and for both alloys. 
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The observed differences between Al and C doping, like the doping range of superconductiv- 
ity, are naturally obtained in the present VGA approach, without the need to invoke other 
factors, as, e.g., interband scattering. These results emphasize that a quantitative descrip- 
tion of the superconducting properties of the two MgB2 alloys require a proper treatment of 
the doping at least on the level of VGA, and suggest that interband scattering plays only a 
minor role. 
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